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We look at multiparticle production processes from the nonextensive point of view. 
Nonextensivity means here the systematic deviations in exponential formulas pro- 
vided by the usual statistical approach for description of some observables like 
transverse momenta or rapidity distributions. We show that they can be accounted 
for my means of single parameter q with \q — 1| being the measure of nonextensiv- 
ity. The whole discussion will be based on the information theoretical approach to 
multiparticle processes proposed by us some time ago. 



1 Introduction 

The multiparticle production processes are usually approached first by means of 
statistical models in order to make quick estimations of such parameters as 
"temperature" T or "chemical potential" fi. The quotation marks reflect the con- 
troversy in what concerns the applicability of such terms in the realm of hadronic 
production. Actually, as has been pointed long time ago [2] , the fact that out of N 
measured particles we usually observe only one of them is enough to make the rele- 
vant distribution (in transverse momentum pt or in rapidity y) being exponential: 



f(X) ~ exp 



X 
A 



(1) 



This is because the remain N — 1 particles act as a " heath bath" to the observed 
one and their action can be therefore summarized by a single parameter A. Iden- 
tification of A with temperature T means that we assume that our system is in 
thermal equilibrium. However, in reality such "heath bath" is not infinite and not 
homogenous (i.e., it can contain domains which differ between themselves). In such 
case single parameter A is not sufficient to summarily describe its influence on the 
particle under consideration. One can regard it as being X-dependent and use 
the first two terms of its expansion around some value Xo (instead of only Xq): 
A ~ Aq + a ■ (X — Xq). In this case distribution (JJJ becomes 



f q (X)~ exp 



X' 


def 


.~Ao. 





l-(l-9) 



X_Xo 
A 



where 



q=l + a. (2) 
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In case when X is energy E and A is temperature T the coefficient a is just the 
reverse heat capacity, a — 1/ Cy 3 . Such phenomenon is called nonextensivity and 
with q being the nonextensivity parameter. Notice that for q — > 1 (or for a — ► 0) 
becomes /(X). 

In next Section we shall look at this problem from the information theory of view 
using Shannon and Tsallis information entropies. Confrontation with experimental 
data will be provided in Section 3. Section 4 will contain summary and conclusions. 

2 Information theory and multiparticle production processes 

Let us consider situation, which happens quite often in physics realm. Suppose 
that experimentalist provided us with some new data. Immediately these data 
are subject of interest to a number of theoretical groups, each of them quickly 
proposing their own distinctive and unique (in what concerns assumptions and 
physical concepts) explanation. Albeit disagreeing on physical concepts they are 
all fairly successful in what concerns fitting data at hand. The natural question 
which arises is: which of the proposed models is the right one? The answer is: to 
some extend - all of theml This is because experimental data are providing only 
limited amount of information and all models are simply able to reproduce it. To 
quantify this reasoning one has to define the notion of information and to do this, 
one has to resort to Information Theory (IT) a . It is based on Shannon information 
entropy 



where Pi denotes probability distribution of interest. Notice that the least possible 
information, corresponding to the equal probability distribution of N states (i.e., 
Pi = l/N), results in maximal entropy, S — hiN. The opposite situation, when 
only one state is relevant, i.e., pk — 1 and Pi^k = 0, results in the minimal entropy, 
5 = 0. Usually one always has some a priori information on experiment, like 
conservation laws and results of measurements, which is represented by a set of 
known numbers, One is then seeking probability distribution {p,} such that: 

• it tells us the truth, the whole truth about our experiment, what means that, 
in addition to be normalized, YliPi = 1> ^ reproduces the known results, i.e., 



• it tells us nothing but the truth about our experiment, i.e., it conveys the least 
information (only the information contain in this experiment). 

To find such {p{\ one has to maximize Shannon entropy under the above con- 
straints (therefore this approach is also known as MaxEnt method) . The resultant 

a The complete list of references concerning IT relevant to our discussion here and also providing 
some necessary background can be found in 14151 . 




(3) 




(4) 



Multiparticle production processes. 



3 



distribution has familiar exponential shape 



1 

Pl = ~z ' exp 



(5) 



Although it looks identical to the " thermal- like" (Boltzmann-Gibbs) formula 
there are no free parameters here because Z is just normalization constant assuring 
that ^2 Pi = 1 and are the corresponding lagrange multipliers to be calculated 
from the constraint equations 

It is worth to mention at this point 4 that the most probable multiplicity 
distribution P(n) in the case when we know only the mean multiplicity (n) and 
that all particles are distinguishable is geometrical distribution P(n) = (n)"/(l + 
(n))( n+1 - ) (which is broad in the sense that its dispersion is D(n) ~ (n)). Additional 
knowledge that all these particles are indistinguishable converts the above P{n) into 
Poissonian form, P(n) = exp(— (n))/nl, which is the narrow one in the sense that 
now its dispersion is D(n) ~ y/ (n). In between is situation in which we know 
that particles are grouped in k equally strongly emitting sources, in which case 
one gets Negative Binomial distribution j7J c P(n) = T(k + n)/ [T(n + l)T(k)] • 

k r .1 k+n 



(<»>) 



1 4- -A- 



/ 

The other noticeably example concerns the use of IT to find the minimal set 
of assumptions needed to explain all multiparticle production data of that time 
[8]. It turned out that all competing models like: multi-Regge, uncorrelated jet, 
thcrmodynamical, hydrodynamical etc., shared (in explicit or implicit manner) , two 
basic assumptions: 

(?) that only part W — K ■ y/s (0 < K < 1) of the initially allowed energy y/s is 
used for production of observed secondaries (located in the center part of the 
phase space); in this way inelasticity K found its justification it turns out 
that K ~ 0.5; 

(ii) that transverse momenta of produced particles are limited and the resulting 
phase space is effectively one-dimensional (dominance of the longitudinal phase 
space). 

All other assumptions specific for a given model turned out to be just spurious^. 



^Notice that using the entropic measure S = . [p;lnpi =p (1 ± pi)ln(l ± Pi)] (which, however, 
has nothing to do with IT) would result instead in Bose-Einstein and Fermi-Dirac formulas: 
Pi = (1/Z) ■ [exp[/3(ei - ft)] =F 1]~\ where (3 and fi are obtained solving two constraint equations 
given, respectively, by energy and number of particles conservation 0. It must be also stressed 
that the final functional form of p; depends also on the functional form of the constraint functions 
Rk{xi). For example, R(x) <x ln(x) and ln(l — x) type constrains lead to pi <x xf(l — Xi)P 
distributions. 

c It is straightforward to check that Shannon entropy decreases from the most broad geometrical 
distribution towards the most narrow Poissonian distribution. 

d The most drastic situation was with the multi-Regge model in which, in addition to the basic 
model assumptions, two purely phenomenological ingredients have been introduced in order to 
get agreement with experiment: (i) energy s was used in the scaled (s/sq) form (with so being 
free parameter, this works the same way as inelasticity) and (ii) the so called "residual function" 
factor e^'* was postulated (t = — (pi — Pj) 2 and (3 being a free parameter) to cut the transverse 
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Suppose now that some new data occur which disagree with the previously 
established form of {pi}. In IT approach it simply signals additional information to 
be accounted for. This can be done either by adding a new constraint (resulting in 
new Afc+i) or by checking whether one should not use exp g (...) rather than exp (...). 
This brings to IT the notion of nonextensivity introduced before and is connected 
with freedom in choosing the form of information entropy. The point is that there 
are systems which experience long range correlations, memory effects, which phase 
space has fractal character or which exhibit some intrinsic dynamical fluctuations 
of the otherwise intensive parameters (making them extensive ones, like T in the 
example above) . Such systems are truly small because the range of changes is of the 
order of their size. In this case the Shannon entropy is no longer good measure 
of information and should be replaced by some other measure. Out of infinitely 
many possibilities we shall choose Tsallis entropy defined as 

^-^Ef 1 -?.)- ( 6 ) 

i 

which, as can be easily shown, goes over to Shannon form J3J for q — > 1. To get the 
{pi} one proceeds in the same way as before but modifying the constraint equation 
in the following way 

i 

This leads to formally the same formula for pi = pf^ as in eq. (J5J but with 
Z — ► Z q and exp(...) — > exp q (...) q . (For our purpose this method is sufficient and 
there is no need to use the more sophisticated approach exploring the so called 
escort probabilities formalism, see . Because such entropy is nonextensive, i.e., 
Sq(A+B) = S q A + SqB + {^ — q)S q A- S q B (see [S] for more details), the whole approach 
takes name nonextensive (Tsallis) statistics. Nonextensivity parameter q accounts 
summarily for the all possible sources of nonextensivity mentioned before, however, 
here we shall be interested here only in intrinsic fluctuations present in the system. 
This is because in this case, as was demonstrated in |12j . parameter q is measure 
of such fluctuations, namely for system described by eq. @ one has that e 

((1/A) 2 ) - (1/A) 2 ]/(1/A) 2 . (8) 

In what follows we shall argue that multiparticle production data on both transverse 
momentum pt and rapidity y distributions can be interpreted as showing effects of 
such fluctuations. 



1 ± 



3 Confrontation with experimental data 

One should stress that what we present here is not a model but only least biased 
and most probable (from the point of view of information theory) distributions 

part of the phase space. Therefore so and /3 were the only important parameters whereas all other 
model parameters were simply irrelevant. 

e Strictly speaking in 1121 it was shown only for fluctuations of 1/A given by gamma distribution. 
However, it was soon after generalized to other form of fluctuations and the word super statistics 
has bee coined to describe this new phenomenon 1131 . 
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describing available data and accounting for constraints emerging either from the 
conservation laws or from some dynamical input (which in this way is subjected to 
confrontation with data). Suppose that we have mass M, which hadronizes into N 
secondaries and hadronization is taking place in 1— dimensional longitudinal phase 
space. Dynamics is hidden here in the fact that decay is into N particles and that 
transverse momenta of particles are limited so one can use their transverse mass 
Ht (suppose that we are not interested for a moment in details of the transverse 
momentum distribution). Once this is settled IT gives us clear prescription what 
to do: 



Look for single particle distribution in longitudinal phase space taken as 
y— space and written as suitable probability distribution, 

, , 1 dN 1 f 2 Edu 
My) = T7 ' -r = Tt / <rPT-n-- (9) 



N dy a N N J ^ d 3 p 
Maximize the following entropy functional 

pY M 

S = - dyf N (y)-ln[f N (y)}, (10) 

J-Ym 

under conditions that /jv(j/) is properly normalized and that one conserves en- 
ergy (because of the symmetry of our problem the momentum will be satisfied 
automatically; no more information will be used): 



f N (y) = l and / dyE- f N (y) = fi t dyf N (y) coshy = — . (11) 

-Y M J-Ym J-Ym 7V 

As result one obtains that the most probable (and least biased) rapidity distribution 
is 

i rY M 

/jv(y) = -7^ ■ exp [— /3 ■ coshy] where Z= I dy exp [— /3 • coshy] (12) 

£ J-Ym 

and where (3 = (3(M, N) is obtained by solving eq. I|ll|) . The y— space is limited to 



ye(-Y M ,Y M ), Y M =ln- 



M'ih 

2^ T 



i + d-MV /2 

' M' 2 J 



, M' = M-(N-2)nr. 



The following points should be remembered: 



• For N = 2 both particles must be at the ends of the phase space, i.e., /a(y) = 
i [8(y - Y M ) + 5{y + Ym)}, therefore formally 0(M, 2) = -oo. 

• Between N — 2 and N ~ No = 2 In N max (where N max = M/fix) P is negative 
and approaches zero only for N = Nq. Therefore only for N — Nq we have 
In(v) = const, i.e., the famous Feynman scaling-^. 



fit was jus a coincidence that at ISR energies this condition was satisfied. But such behaviour 
of N as function of energy was only transient phenomenon and there will never be scaling of this 
type, notwithstanding all opinions to the contrary heard from time to time. 
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• For N > Nq additional particles have to be located near the middle of phase 
space to in order to minimize energy cost of their production. As result (3 > 
now, in fact (see @] for details) for some ranges of < E >= M/N quantity 
(3 — (3 < E > remains approximately constant. Needles to say that for N — ► 
N max all particles have to stay as much as possible at the center and therefore 
(3(M,N ->■ N max ) -> +oo whereas fN->N ma Av) ~> %)■ 





y -Is [GeV] 

Figure 1. Left panel: examples of rapidity distributions for different nonextensivities. Right panel: 
the energy dependence of the NB parameter 1/k (line) together with results for the nonextensivity 
| q — 1 1 obtained when fitting pp and pp data at different energies (points) - see the top-left panel 
of Fig. 



For nonextensive approach the only changes to be performed is to replace exp(...) 
by expj...) q in the way explained before. In addition N — > Nq^ = (2 \n.N max ) q . 
The changes in the shapes of f q (y) function for different values of q are shown 
in Fig. . Notice that for q > 1 one enhances tails of distribution and at the 
same time reduced its height. For ql the effect is opposite (in this case there is 
kinematical constraint on the allowed phase space, namely it must conform to the 
condition that 1 — (1 — q)j3 q ^T coshy > 0. In fact in JS] we have found that we 
can fit reasonably well data on pp |16j with q being the only parameter. Resultant 
q < 1, which was cutting off the available phase space was there playing effectively 
the role of inelasticity of reaction mentioned before' 1 . 

We repeated attempt to fit pp (and this time also pp) data JHI again but this 
time accountig properly for inelasticity (which was now free parameter to be de- 
duced from data) and also for other possible fluctuations which hadronizing process 
offers us and which were represented by parameter q. We shall not discuss here the 
inelasticity issue (see ^T] for details and further references) but concentrate on the 
q parameter. In Fig. 2] we show the values of this parameter obtained for different 
energies compared with the values of the parameter k of Negative Binomial dis- 
tribution (NB) taken for different energies from existing data ^7]. The agreement 
is perfect what prompted us to argue that parameter q is accounting for a new 

9 In general situation is bit more involved, see ,for example, 1141 for more details. This, however, 
is out of scope of our presentation. 

h The aim of 1151 was to provide the cosmic ray physicist community a kind of justification of 
the empirical formula they used, namely that f(x) oc (1 — a ■ x) n , where x is Feynman variable, 
x = 2E/^/s, and where a and n are free parameters. They turned to be both given by q only 1151 . 
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Figure 2. The most characteristic examples of comparison with data on rapidity distributions for 
different reactions. Top-left panel: the best fits to pp and pp data | 16 | using nonextensive approach 
(cf. Illl V The resulting values of the q parameters are shown above in Fig. Top-right panel: 
fits to NA49 data for negatively charged pions produced in central PbPb collisions at different 
energies 1181 (cf. The best fit for 17.3 GeV is, however, for two extensive sources located 

at y = ±Ay = 0.83 in rapidity (total inelasticity in this case is K = 0.58). Bottom-left panel: 
fits to PHOBOS data for the most central Au + Au collisions 1191 (cf., 1201 ), Bottom-right panel: 
rapidity spectra measured in e+e — annihilations at 91.2 GeV 1211 (dotted line is for K q = 1 and 
q = 1 whereas full line is our fit with K q = 1 and q = 0.6, see 1201221 1. 



8 



O.V.Utyuzh, G.Wilk, Z.Wlodarczyk 



bit of information, so far unaccounted for, namely that in addition to the mean 
multiplicity which was our input parameter there are other data |17j showing the 
whole multiplicity distributions, in particular showing their NB character. Actually 
this conjecture is supported by the observation |7] that NB can be obtained from 
Poisson distribution once one allows for fluctuations in its mean value n proceeding 
according to gamma distribution, namely 

_ [°° j_ e-*n n 7 fc n fc - 1 e-^ _ Tjk + n) 7 " n ^ 

[n) ~ J n\ ■ r(fc) r(i + n)r(fc) ' ( 7 + i) fc +» 1 ' 

where 7 = -Mv and \ = a , ^ n °$ — , 1 , . Assuming now that these fluctuations 
contribute to nonextensivity defined by the parameter q, i.e., that D(n) = q — 1 
one gets that 

= 1 + J, (14) 

what we do observe (cf. also \'2'6\). 

The next two panels show comparison with AA data taking at SPS (NA49) and 
RHIC (PHOBOS) energies. The NA49 data can be fitted with q = 1.2, 1.16 and 
1.04 going from top to bottom and so far we do not have plausible explanation for 
these number. However, it should be stressed that at highest energy two component 
extensive source is preferable. The PHOBOS data with q = 1.27, 1.26 and 1.29 
going from top to bottom. Extensive fits do not work here at all. The example 
of e + e _ data clearly show that in this case there is some new information we did 
not accounted for. In this case the whole energy must be used but only q < 1 fit 
(cutting-off a part of phase space) can be regarded as a fair one. But even then we 
cannot obtain minimum at y = 0. It looks like we have two sources here separated 
in rapidity (two jets of QCD analysis) but of no statistical origin (rather connected 
with cascading) |22| . 

Let us now proceed to examples of pr distributions, see Fig. [3J Both pp data at 
high energies and nuclear data at lower energies are shown. They all can be fitted 
only with q > 1 (for details see and 121]. The characteristic feature is that now 
values of q are much smaller than those obtained fitting data in longitudinal phase 
space. We shall not pursuit this problem here referring to discussion in [5| . Instead 
we shall concentrate on the nuclear data (right panel) and stress here that these 
data can, in our opinion, be connected with fluctuations of temperature mentioned 
at the beginning. In fact the q = 1.015 obtained here fluctuation of T of the order 
AT/T ~ 0.12, which do not vanish with increasing multiplicity In fact they are 
fluctuations existing in small parts of the hadronic system with respect to the whole 
system rather than of the event-by-event type, for which AT/T ~ 0.06/ViV — > 
for large N. Such fluctuations are potentially very interesting because they provide 
a direct measure of the total heat capacity of the system, C, because 

* 2 W) def 1 

-W = c = q - 1 - (15) 

In fact, because C grows with the volume V of reaction we expect that 
q(hadronic) > q(nuclear) which seems to be indeed observed (cf., for example, 
Fig. EJ. 



Multiparticle production processes. 



9 




Figure 3. Transverse momentum spectra showing nonextensive character. Left panel: fits to pt 
spectra from UA1 experiment on pp collisions [241 . Here q = 1.11, 1.105 and 1.095 for energies 
decreasing from 900 to 200 GeV, respectively (cf., 1111 ), Right panel: fit to pj< spectra from S + S 
collisions 1251 (cf., 1201 '). Notice that in this case q = 1.015 only, i.e., it is smaller less than for pp 
case, see eq. 1151 . 



4 Summary and conclusions 

We have demonstrated that large amount of data on multiparticle production can be 
quite adequately described by using tools from information theory, especially when 
allowing for its nonextensive realization based on Tsallis entropy. We have argued 
that nonextensivity parameter q entering here can, in addition to the temperature 
parameter T of the usual statistical approaches, provide us valuable information on 
dynamical fluctuations present in the hadronizing systems. Such information can 
be very useful when searching for phase transition phenomena, which should be 
accompanied by some special fluctuations of nonmonotonic character. It is therefore 
worth to end with saying that some interesting fluctuations of this sort seems to be 
already seen and investigated [2Z|. All that calls for some more systematic effort 
to describe existing data in terms of (T, q) for different configurations and energies 
in order to find possible regularities in their system and energy dependencies and 
possible correlations between them. 
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